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A  NOTE  ON  THE  PROPAGATION  OF  CERTAIN  LF  PULSES 
UTILIZED  IN  A  RADIO  NAVIGATION  SYSTEM 

by 

J.   R.   Johler 

ABSTRACT 

Low  frequency  pulses  utilized  in  a  radio  navigation  system  can 
be  synthesized  and  analyzed  by  a  specialization  of  the  general  theory 
of  pulse  propagation.     The  analysis  technique  required  to  accomplish 
such  a  research  is  detailed,    and  the  introduction  of  receiver  and 
detector  circuits  into  such  an  analysis  is  described.     The  time -dif- 
ference dispersion  correction,    or  the  "discrepancy  correction"  for 
operation  of  a  radio  navigation  pulse  system  on  either  the  ground- 
wave  time  mode  or  the  ionospheric  time  modes  of  propagation  is 
defined. 


A  NOTE  ON  THE  PROPAGATION  OF  CERTAIN  LF  PULSES 
UTILIZED  IN  A  RADIO  NAVIGATION  SYSTEM 

by 

J.    R.    Johler 

The  interest  in  recent  years  in  low  frequencies  (LF)  has  been 
inspired  by  the  development  of  the  Loran-C  (Cytac)  [Frantz,    Dean, 
and  Frank,    1957  ]  and  other  [Schneider,    1955]  [Ross,    1959]radio 
navigation  and/or  timing  systems.     The  high  reliability  and  pre- 
dictability requirements  of  such  systems  leads  the  researcher  and 
system  designer  quite  naturally  to  low  frequencies. 


L1NE_0F  — POSITION  ,    At    =    CONSTANT 


Figure  1.     Radio  navigator1  s  hyperbolic  line -of -position. 


The  principle  of  radio  navigation  can  be  described  by  the  line -of  - 
position,    figure  1.     A  master  transmitter,    M,    and  a  slave  "transponder", 

S,    are  separated  by  the  length  of  baseline,    d.  ,    which  is  ordinarily  a 

b 

constant.     An  observer,    O,    with  a  phase -difference  indicator  com- 
pares the  phase  of  a  signal  transmitted  to  O  directly  (i.  e. ,    a  trans- 
mission with  which  the  distance,    d     ,    is  associated)  with  the  signal 

m 
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received  via  the  slave,    S,    (i.  e, ,    a  transmission  with  which  the 

distance,    d,  ,    the  slave,    S,    and  the  distance,    d   ,    are  associated), 
b  s 

If  the  phase-difference,    Acj>,    can  be  interpreted  as  a  time,    in  micro- 
seconds,    (t  =  10     cj>/u>  at  a  transmission  frequency,    f  =  oo/Zir)  then  the 
time -difference,    At,    represents  the  phase  measurement,    or  phase 
lock,    of  the  automatic  system. 


At  =      ^1     [d,    +  d     -  d     ]  +  C     ,  (1) 

— —  L    b  s  rn  s 


(where  C     is  the  slave  (S)  coding  delay  which  is  a  system  constant, 

and  r)  /c  is  the  reciprocal  of  the  speed  of  light).     The  time  difference, 

At,    as  a  constant,    together  with  a  constant  baseline,    d   ,    describes  a 

hyperbolic  locus  of  points  known  as  a  line -of -position,   figure  1.     Two 

such  lines -of -position  (which  can  be  generated  with  the  aid  of  one 

additional  slave)  determine  a  navigation  "fix".     The  transmission 

over  the  distances,    d,  ,    d   ,    d     ,    is  accomplished  by  use  of  LiF  radio 

b        s        m 

wave  propagation.     It  is  therefore  evident  that  the  system  designer 
is  interested  in  the  L.F  field,    E(oj,  d)  volts  /meter,    and  in  particular 
the  phase  of  the  field,    arg  E(oo,  d),    at  frequencies,    f  =  w/2tt  propagated 
over   these  distances.     The  existent  Loran-C  system  employs  the 
ground-wave  mode  [Johler,    Kellar,    and  Walters,    1956],    [Johler, 
Lilley,   and  Walters,    1959],    [Johler,    Lilley,    and  Walters,    19  60]  for 
transmissions  over  these  distances  to  the  exclusion  of  the  ionosphere 
modes.     Special  instrumentation  which  utilizes  pulse  techniques  and 
special  radiated   pulse  signals,  E(t'  ,  d)  volts/meter,    which  can  be 
readily  tagged  at  a  point-in-time  where  t'    is  reckoned  locally  from 
the  beginning  of  the  pulse,    t1    =  t  -  r\  d/c,    is  required  to  achieve  the 
ground-wave  "phase  lock"  to  the  exclusion  of  the  ionospheric  modes. 
Thus  it  is  necessary  to  consider  the  influence  of  the  conduction  and 


-3- 

displacement  currents  in  the  ground,    the  curvature  of  the  ground, 

and  the  proximity  of  the  source  (induction  and  electrostatic  fields). 

It  is,    therefore,    necessary  to  develop  and  apply  a  time  correction, 

t     =  t   (d),    for  each  distance  transmission  associated  with  the  system, 
c        c 

The  time -difference  equation  is  therefore  corrected, 


At  =      ^1      [±    +  d     -  d     ]  +  t   (d,  )  +  t   (d   )  -  t   (d     )    +  C  (2) 

- —  L    b         s         mJ         c     b  c     s  cm  s    , 


to  describe  the  line -of -position  with  greater  accuracy. 

Recent  proposals  [  Doherty,    Hefley,    and  Linfield,    I960]  have 
suggested  the  employment  of  the  world-wide  Loran-C  system  to  move 
time  around  the  world  by  means  of  such  a  "phase  lock"  to  an  accuracy 
of  one  microsecond  or  better.      The  ability  of  the  system  to  synchronize 
to  fractional  microsecond  precision,    and  the  corresponding  predicta- 
bility of  the  ground-wave  mode  of  propagation  imply  the  feasibility 
of  this  application,    hence  the  notation:      Loran-C  radio  navigation 
timing  system. 

Since  the  Loran-C  radio  navigation  system  employs  pulse 
techniques,    a  discussion  of  the  analysis  procedures  requires  special 
consideration  of  the  pulse  type  transmission,    E(t'  ,  d)  volts /meter. 
The  general  approach  [  Johler  and  Walters,    1959]  was  developed  so 
that  the  CW  field,    E(oo,  d),    could  be  related  to  the  transient  field, 

E(t'  d),    and  the  source,    F    (t),    by  a  direct  evaluation  of  the  Fourier 

s 

transform  integral, 

00  00 


E  (t«,  d)  =  ~      {     exp  (ia)t»)  E  (W|  d)     f    exp  (-iu>t  )  Ffl  (t)  dt  do 


(3) 


-00  w 

00 

or      E  (tSd)  =  ~      I    exp  (icjt»)  E  (w,  d)  Fg  (u>)  doj    .  (4) 

-00 
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The  complex  source  function,    F    (t),    introduced  to  simplify  the 
operational  techniques  of  analysis,    also  provided  a  description  of  an 
"ideal"  amplitude  envelope  for  the  pulse,     |E(t'  ,  d)|,    which  was  propa- 
gated simultaneously  with  the  pulse. 

The  transform  of  the  propagation  medium,    E(u),  d)  volts /meter, 
is   a  very  complicated  function  which  has  been  extensively  described 
in  the  literature  [Bremmer,    1949].     The  transform  is  frequently 
represented  geometric -optically  at  LF,    such  that  the  total  field, 
E(oj,  d),    for  the  particular  case  of  a  vertically  polarized  dipole  source 
can  be  written  as  the  sum  of  time  modes  [  Johler,    1961]  of  propaga- 
tion,  j  =  0,    1,    2,    3  •  •  • ,    between  the  ionosphere  and  the  earth, 

P 
E(w,  d)=     Y     E  (o>,d)      ,  (5) 

j=0 

in  which  the  zero  order  (j  =  0)  time  mode,    E    (go,  d),    is  the  ground 

o 

wave.     Each  time  mode  of  propagation,    j  =  1,    2,    3  •  •  •    or  j  ^=  0,    is, 


E.(w,  d)  =  iwd  D."1  C  exp  (iutc.)  G1.  Gr.  a.  F.  C.      ,  (6) 

j  J  J         J        J      J      J      J 


where, 


C  =  Ioib2/4TTKd3  =  (io"7)/d,    (IQi=    1)  (7) 


and  the  local  time  for  the  ground  wave,    t1     ,    is, 

o 


t*     =  t-b      ,  (8) 

o 

b  =  ti^/c     ,  ^j 
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where  [Johler,    Kellar,    and  Walters,    1956] 


E   (co,  d)  =  JE    (co,  d)  j  exp  [  -i(a)b  -1-  $   )]      ,  (10) 

o  o  c 

Q 

where  c  is  the  speed  of  light,    c  ~  3(10    )  meters /second,    and  r\    is 
the  index  of  refraction  of  air,    T)    ~  1.  000338  typically,    and    similarly, 
the  local  sky-wave  time,    or  sky-wave  delay,    t'  .  (j  =  1,    2,    3  •  •  •  )  is, 

t*.=  t-b.      ,  W 

J  J 

b    =  ti    D/c,    (j  a  1,    2,    3««     )  (12) 

J  J 


The  quantity,    b.  -  b,    is  frequently  called  the  relative  sky-wave  delay, 

figure  2  (relative  to  the  ground  wave).     The  physical  length  of  the  ray, 

D.,    can  be  evaluated  geometrically  for  a  reflection  at  an  altitude,    h, 

J 
above  the  surface  of  the  earth  of  radius,    a, 


D.  =  2j  T(a+h)  cos  <j>.    .  -  a  cos  t.]      , 


(13) 


where  cj>.  is  the  angle  of  incidence  of  the  "ray"  on  the  ionosphere  and 
T.  is  the  corresponding  angle  of  incidence  on  the  Earth  (the  subscript 
"j"  reminds  the  reader  that  the  equation  refers  to  the  particular  time 
mode  under  consideration).      The  angles,    <$>.  and  T.,    are  evaluated 
quite  simply  from  the  geometry  if 


•     x  a   -1         •        9  (14) 

sin  <b.    .  =  A.       a  sin    -^-r-       ,  K^' 

cos  <|>        =  A  -1  [a(l-cos  9/2j)  +  h]      ,  (15) 


sin  t.  =  A."    (a  +  h)  sin  0/2j      , 

cos  t.  =  A."    [a(cos  8/2j  -  1)  +  h  cos  9/2j)      , 
J  J 

1 

A.  =  [2a(a  +  h)  (1  -  cos  e/2j)  +  h2]  , 


(16) 
(17) 
(18) 


where  8  is  the  angle  at  the  center  of  the  earth  subtended  by  the  dis- 
tance,   d,    along  the  surface  of  a  spherical  earth,    or,    simply,    d  =  a0. 

t  r 

The  factors,    G  .  and  G  .,  refer  to  the  transmitter  and  receiver 

J  J 

vertical  antenna  patterns  respectively,    and  again  the  subscript,    "j", 

is  employed  to  designate  the  particular  time  mode  under  consideration. 

t  r 

The  vertically  polarized  Hertzian  dipoles  reduce  to  G  .  =  G    .  =  sin  T.. 

Of  course  other  antenna  patterns  can  be  introduced  by  these  factors. 

The  factor,    F.   [Johler,    1961b],    [Wait,    1958],    describes  the 

antenna  pattern  corrections  of  the  transmitter  and  receiver  antennas 

at  grazing  angles  of  incidence  on  the  earth  and  the  factor,    C., 

describes  the  effective  reflection  coefficient  for  multiple  hops  or  a 

single  hop  between  the  ionosphere  and  the  Earth  [Johler,    1961b], 


C,  = 


1 


J'jTX 


cbr1 


1  +  A     x 

I 

1  —  J\  —    X    —  J\  ^    X 


(19) 


x  =   0 


At  =   -R       T 
1  m      mm 


A-  =  R     T        -»-  R       T 
Z         e      ee  m      mm 


A,  =  R     R 
i         em 


-T       T  +  T         T 

ee      mm  em      me 


(20) 


where  R     and  R       are  the  ground-reflection  coefficients  for  vertical 
em 
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and  horizontal  polarization  and  T      ,    T        ,    T        ,    and  T         ,'are  the 
c  ee        em        me  mm 

ionosphere  reflection  coefficients  [  Johler  and  Walters,    I960]  [  Johler, 
1961c]  [Johler  and  Harper,    1962]. 

Since  plane -reflection  coefficients,    T,    are  employed  in  this  type 
of  analysis,    the  convergence-divergence  factor,  <x,   is  necessary  to 
convert  plane-reflection  coefficients  to  spherical-reflection  coefficients, 
or  from  the  view  point  of  ray  theory  it  is  necessary  to  take  account 
of  ray  focusing  by  the  ionosphere  and  a  corresponding  defocusing  of 
rays  by  the  Earth. 

As  a  consequence  of  the  time  modes, 

p  p  00 

e  (t»,  d>  =  y  ej  Wy  d)  =  2  -h  j  exp  ^y  ej  (a)» d>  Fs  (w)  **• 

j=o  j=o         -oc  (21) 

This  merely  represents  the  sum  of  separate  Fourier  integrals  for 
each  time  mode,    separated  in  time  by  the  sky-wave  time -mode 
delay,    t1  .,    figure  2. 
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Figure  3.       Observations  of  time  modes  with  the  aid  of  pulses  of 

various  lengths  illustrating  the  separation  of  time  modes 
and  the  measurement  of  a  point-in-time  on  the  leading 
edge  of  the  pxilse  with  the  "amplitude  envelope  minus  the 
derivative  of  the  amplitude  envelope"  method,    (August 
1953,  f  =  100  kc/s). 

The  sum  of  the  time  modes  is  illustrated     for  various  (Cytac) 
Loran-C,    Radio  Navigation  System  pulses,    figure  3.     The  undisturbed 
form  of  the  radiated  pulse  is  shown  by  oscillograms  which  were  made 
at  short  distances.     The  action  of  the  propagation  medium  is  evidenced 
by  the  complicated,    multiple,    superposed  pulses  observed  at  great 
distances.     Thus  the  oscillogram  or  amplitude -time  function, 
Re  E(t'  ,  d),    observed  at  great  distance,    comprises  the  superposed 


From  the  author'  s  personal  notes  on  the  1953-1954  preliminary 
testing  of  the  Loran-C  (Cytac)  system. 
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time  modes,    E.(t'  ,  d),    j  =  0,    1,    2,    3  '  '  ' ,    apparently  added  together 

according  to  the  complex  field  specification,    E.,    of  each  and  separated 

in  time  by  the  sky-wave  delay,    t1  .,    of  each  time  mode,   j.     The 

approximate  area  in  the  amplitude -time  plane  occupied  by  each  mode 

of  propagation  is  shown,    figure  3.     Areas  of  overlap  in  which  both 

constructive  and  destructive  interference  between  cycles  of  the  multiple 

pulses  produces  notches  and  enhancements  in  the  composite  pulse 

can  be  observed. 

The  instrumentation  which  has  been  developed  for  the  Cytac  and 

Loran-C  systems  is  quite  capable  of  elegantly  tagging  cycles  and 

minute  fractions  of  a  cycle  on  each  of  the  multiple  pulses.     The 

envelope  detector,    in  effect,    forms  the  difference  between  the  amplitude 

envelope,     |E(t',d)|,   and  its  derivate,  ;     JE(t',d)|,   and  an  elaborate 

servomechanism  seeks  out  a  null  point  or  zero  crossing  time  for  the 

time -function,    thus  determining  a  time -root,    T    ,    of  the  differential 

c 

equation,    illustrated,    figure  3, 

F(t',d)=c1    |E(t\d)|  -  c2-lrT     |E(f,d)|=0,  (22) 

where  c    and  c     are  constants  which  move  the  point-in-time  which 
is  to  be  tagged  on  the  oscillogram  and  which  can  be  set  by  the  operator 
to  pick  out  a  particular  point  on  the  pulse.       Both  amplitude  and  time 
(phase)  can  be  measured  at  such  a  point.     Additional  techniques  to 
detect  cycles  and  minute  fractions  of  a  cycle  under  the  envelope  on 
the  leading  edge  of  a  pulse  have  also  been  developed  [Frantz,    Dean, 
and  Frank,    1957],    and  the  appropriate  time-correction,    t   ,    of  the 
continuous  wave,    E(co,  d),    is  appropriate  in  this  case. 
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The  earliest  time-root,    T     =  T     (d),    of  this  equation,    F(t'  d)  =  0, 

c  c 

constitutes  the  dispersion  correction  for  the  navigator1  s  line-of- 
position  (1,    2)  for  propagation  via  ionospheric  modes,    or 

At  =  l±      [D.(db)+D.(ds)-D.,dm)]+Tc<db)  +  Tc(ds) 

-Tc(dm1  +  CS      .  <"> 


The  ground-wave  mode  is  characterized  by  negligible  dispersion 
for  the  type  of  pulses  employed,    T    (d)  =  t    (d)  +  c,    where  c  is  a  con- 
stant which  becomes  part  of  the  coding  delay,    C    ,    and  of  course,    the 
ray  length,    D.(d)  =  d,    by  definition. 

The  theory  which  can  predict  the  behavior  of  pulses  is,    therefore, 
an  important  adjunct  to  a  system-prediction  analysis.     The  radiated 
field  from  a  step  or  impulse -cur  rent  dipole  source  has  been  the  sub- 
ject of  considerable  discussion  among  theoreticians   [  Johler  and 
Walters,    1959]  since  such  functions  are  quite  simple  to  evaluate. 
Thus  the  field  of  a  step-function  dipole  current  source  can  be  written 

CO 

E(t*,  d)  =  i      f    ]E(W|  d)  ]  J  i  sin  [art*  -  tfj  j    dco      ,  (24) 

o 

TT 

and  the  particular  case  of  the  ground-wave  mode,    writing  (j)1     =  <j>     -  — 

E    (go,  d)  =    |E    (oj,  d)|   exp  [  -i  <j>    1    is  illustrated,    figure  4,    (d  =  1000 
o  o  L  c 

statute  miles).      The  impulse-source-current  dipole  (vertical  polariza- 
tion) reduces  to  the  following  simple  formula, 

E(t*,d)  =  ^     ^    ]E(w,d)  |  jcos  [cot*  -♦»]}    du      .  (25) 
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The  results  of  a  quadrature  evaluation  [  Johler  and  Walters,    1959] 
[  Johler  and  Lilley,    1961]  of  this  integral  are  also  illustrated,    figure 
4,    for  the  ground-wave  mode  of  propagation. 


t 1 1 r 


VERTICAL  ELECTRIC  DIPOLE  SOURCE 
Ioi  =  l  ampere  meter  (step  function) 
d  =  1000  miles 
cr  =  0.005  mhos /meter 
£,  =  15 


Figure  4.     The  step  and  impulse -source  current  radiation 
at  a  great  distance  from  the  source. 


These  pulses  give  physical  insight  into  the  propagation  mechanism, 
but  are  of  little  assistance  to  the  system  designer.     More  sophisticated 
pulses  must  be  developed  mathematically  and  the  transfer  function  of 
the  system  equipment  must  be  considered. 

The  finite  bandwidth  of  a  receiver  or  filter  can  be  introduced, 
equations  (3  and  4),    by  a  mutilation  of  the  Fourier  transform  with 
one  or  more  frequency-selective  networks.     The  intervals  can  then  be 
evaluated  by  quadrature  methods  [Johler  and  Walters,    1959].     In 
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principle,    this  simplifies  the  numerical  problem  of  integration,    since 
the  range  of  integration  necessary  to  achieve  convergence  of  the 
integral,    already  limited  by  the  propagation  medium,    is  further 
limited  by  the  operation  of  the  receiver  transfer  function  or  filter  on 
the  pulse  image.     The  generalized  Fourier  transform,    for  some  com- 
plicated source  function,    can  then  be  represented  as  follows: 

f(co,d)=  E(co,d)  [fSsl^)  +  fSj>2^)  +  0"][^8l(w)fr82^)---].  (26) 

f      ,(03)  f        (oj)  •  •  •  are  source  transforms  which  have  been  superposed 
s,  1  s,  c 

to  determine  some  complicated  waveform,     f      ,(<jj),    f        (co)  •••   are 

r,  1  r,  d 

transfer  functions  of  the  frequency-selective  network  which  have 

been  introduced  to  mutilate  the  transform  or,    in  effect,    restrict  the 

frequency  spectrum  of  the  pulse. 

The  damped  sinusoid,    F    (t)  =  exp  (-vt),    v  =  c,  +  ioo  ,    where  c. 

s  1  c  1 

is  the  damping  and  f     =  cj   /2tt  is  the  characteristic  frequency,    is  a 
more  sophisticated  pulse  [  Johler  and  Walters,    1959]  which  has  found 
useful  application  to  system  studies.     It  is  necessary,    however,    to 
take  into  consideration  the  receiver  or  filter.     Perhaps  the  simplest 
case  is  the  following  series  network: 


o 1| r^nnp— 

C  L 


R 


The  frequency  response,    f    (s),    s  =  iw,    of  this  network  can  be  written: 

(A)  s 

8       + 8    + 


m 


CL 
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The  damped  cosine  source  model  [Johler  and  Walters,    1959], 
exp  (-vt),    v  =    c    +  iw     ,    c 
form,    f(s),    is  as  follows: 


exp  (-vt),    v  =    c    +  iaj     ,    can  be  assumed  in  such  a  way  that  the  trans 


*(s) 


(<*  4-  P)    a 
(s+vM8+a)(s+p) 


(28) 


or, 


f(s)  = 


off 


A(s+v) 


+ 


v  a 


-  i 


A(s+a  ) 


A  (s+p) 


where  A  is  described  by  the  determinant, 
1  1  1 


A    = 


and, 


a  - 


v+p 


ap 


2L 


v+a 


a+p       a+p 


va 
ap 


(29) 


R  v 
L 


CL 


(30) 


P    = 


2L 


CL 


(31) 


For  the  special  case, 
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^2< 


L  '  CL 


j3  =  a  *  (32) 

i.  e.  ,   a,  and  (3  are  conjugate  frequencies.     This  special  case  character- 
izes a  "ringing"  circuit  as  a  result  of  the  introduction  of  the  oscillatory 
terms  with  frequencies,    Im  a  and  Im  (3  ,   which  correspond  to  the  resonant 
frequency  of  the  circuit.     Exponential  damping,    exp  (-  R/2Ij),    is  also 
introduced  by  the  circuit.     It  can,    therefore,    be  concluded  that  the 
effect  of  a  simple  series -resonant  circuit  is  the  same  as  the  sum  of 
three  oscillatory  waves  as  follows: 


E(t«,d)    =    A    Ev(t»,d)  +  B    Ea(t»,d)  +  C  E3  (t»,d) 


(33) 


where  each  wave  is  calculated  by  published  [  Johler  and  Walters, 
1959]  quadrature  methods  at  the  following  complex  frequencies: 

v    =  c,   +  iw 

1  c 

a   =  c,  +  iw 

2  r 

P    =  C3  "  itJr  (34) 


where, 
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c2=  c3  = 


2L 


R_ 
L 


CL 


c3    =    3 


-      A  -*     \    T. 


CL 


U) 


R_  \2        _4_ 
L  <*  CL 


ft) 


R^2<  JL 


CL 


w      =  0 

r 


ft 


2,  ± 

CL 


(35) 


The  constants  A,    B,    and  C  are  defined  as  follows! 


A=A 


-1 


v(3  VCL 

a(3         af(3 


B  =A 


■if"  vQ 


C  =A 


-1 


1  - 


v(3 

a]3 


(36) 
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Figure  5.      Comparison  of  the  calculation  of  the  pulse 

observed  through  an  infinite  bandwidth  receiver 
with  the  calculation  of  the  same  pulse  observed 
through  a  finite  bandwidth  receiver,    illustrating 
the  effect  of  'narrow  band"  as  compared  with 
"broad  band"  receiver. 

The  more  complicated  filters,    f   (<*>),    are  introduced  as  a 

r 

mutilation  of  the  Fourier  transform,    f(w,  d),    simply  by  complex 
multiplication, 


f(u>,d)=  E(w,d)fs(co)fr(u>) 


(37) 


The  filter  transfer  function,    f   (w),    can  be  measured  experimentally, 

r 

preferably  at  precalculated  Gaussian  frequencies,   co      [Johlerand 
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Walters,    1959],    suitable  for  quadrature  methods  provided  that  both 
amplitude  and  phase  are  measured  for  each  frequency.     The  Fourier 
transform  (37)  can  then  be  formed  and  the  Fourier  integral  (3,    4)  can 
then  be  evaluated.     The  results  of  such  a  calculation  are  illustrated, 
figures  5  and  6,   for  two  damped  sinusoids  by  a  comparison  of  the 
pulse  in  an  "infinite"  bandwidth  receiver  with  a  receiver  of  restricted 
bandwidth  and,   as  might  be  expected,   the  ringing  of  the  filter  as  a 
result  of  the  oscillatory  terms  (33)  is  evident. 


VERTICAL   ELECTRIC  DIPOLE   SOURCE 

fc  =  100  kc 
c,*  2.5(10*) 
v  =  0.005  mhos  /meter 
d=50  miles 


Time,  t,'  Microseconds 


Figure  6.     Comparison  of  the  calculation  of  the  pulse 

observed  through  an  infinite -bandwidth  receiver 
with  the  calculation  of  the  same  pulse  observed 
through  a  finite -bandwidth  receiver. 


A  useful  model  for  the  type  of  pulse  normally  employed  in  radio 
navigation  systems  can  be  readily  synthesized  by  superposition  of 
damped  sinusoids.     For  example,    the  "classical"  sine-squared  pulse 
can  be  evaluated  by  redefining  the  source  function, 


-18- 


7 
Re  F   (t)  =  exp  (-c.t)  sin    u>  t    sin  w  t,  (38) 

Si  D  c 


or, 

l  1 

Re  Fg(t)  =    j  exp  (-c^t)  sin  u>ct  -  %■  exp  (-Cjt)  sin  (a>c  +  2oop)t 

-  4  exp  (-c.t)  sin  (co     -  2  oo    )t      , 

4       P         1  c  P  (39) 

or,    the  complex  source  function, 

Fg(t)  =    \  exp  (-vt)  -  \  exp  (-v^)  -  i  exp  {-vfl      ,  (40) 

can  be  employed  for  envelope  synthesis,   where, 
v  =  c.  +  ico     (as  before) 

v    =  c.  +  i  (eo     +  2  oo    ) 
lie  p 

v,=  c.+i(u>-2u>),  (41) 

2         1  c  p 

Typically,    the  radio  navigation  type  pulse  has  a  value  for  c    of  the 
order  of  oo  .     It  is  quite  obvious,    however,    that  the  transient  field, 
E(t'  ,  d),    for  the  damped  sine-squared  source-current  pulse  is  merely 
the  sum  of  three  waves  calculated  as  previously  described  damped 
sinusoids, 


E(t'd)  =  -^Ev(t«,d)  +^-  Ev(t«,d)  +  ■£  Ev   (t*,d)      . 

1  C* 


(42) 
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Figure  7.        The  damped  sine-squared  pulse,   illustrating  a  radio 
navigation  system  type  of  pulse  with  corresponding 
amplitude  envelope. 
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Figure  8.       The  damped  sine -squared  pulse,   illustrating  a  radio 
navigation  system  type  of  pulse  with  corresponding 
amplitude  envelope. 
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The  results  of  a  radio  navigation  system  pulse  synthesis  are 
illustrated  for  propagation  of  a  damped  sinusoidal  sine-squared 
source  current  to  a  distance  of  50  to  1000  miles,    figs.    7,    8.     Both  the 
"oscillogram"  or  observed  waveform,    Re  E(t' ,  d),    and  the  synthetic 
amplitude  envelope,     |E(t' ,d)|,    are  illustrated. 
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